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On the basis of the perturbative QCD and the relativistic quark model we investi- 
gate the relativistic and bound state effects in the production processes of a pair of 
5-wave doubly heavy mesons with opposite charge conjugation consisting of b and 
c quarks. All possible relativistic corrections in the production amplitude including 
the terms connected with the transformation law of the bound state wave func- 
tion to the reference frame of the moving pseudoscalar V— and vector V— mesons 
are taken into account. We obtain a growth of the cross section for the reaction 
e + + er — > J '/ty + 7/ c due to considered effects by a factor 2 -j- 2.5 in the range of the 
center-of-mass energy yfs = 6 -r- 12 GeV. 
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I. INTRODUCTION 

Different processes of the production and decay of heavy mesons consisting of heavy b 
and c quarks provide the means for revealing the role of the color and spin quark forces. 
The aim of many present experiments consists in the increase of experimental accuracy what 
is important for the detailed comparison of different existing theoretical approaches to the 
heavy quark problems 0, 0, 0, 0| • The exclusive production of a pair of doubly heavy mesons 
with c-quarks in e + e~ annihilation has attracted considerable attention in the last years. 
This is due to the fact that the cross section of the process e + + e~ — > J/^/ + r\ c which was 
measured in the experiments on Babar and Belle detectors at the energy \fs = 10.6 GeV 

s „, , ,n f 25.6 ±2.8 ±3.4 [5,61 

a(e+e" - J/^> + Vc ) x B(rj c -> 2 charged) = | 1? 6 ± 2 8 +i.5 [? j (1) 
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leads to a discrepancy with the theoretical calculation in the framework of nonrelativistic 
QCD (NRQCD) by an order of magnitude [l, El 0, HI & This conclusion is based on 
calculations in which the relative momenta of heavy quarks and bound state effects in the 
production amplitude were not taken into account. A set of calculations was performed 
to improve the nonrelativistic approximation for the process. In particular, relativistic 
corrections to the cross section a(e + e~ —>■ J/^>n r ) were considered in a color singlet model 
in Ref.pl using the methods of NRQCD It was obtained here that the relativistic 

corrections increase the cross section by a factor 2.4 for the production J/^> + r\ c . Another 
attempt to take into account the relativistic corrections was done in the framework of the 
light-cone formalism Q Q Q. Here it was shown that the discrepancy between the 



experiment and the theory can be eliminated completely by considering the intrinsic motion 
of heavy quarks forming the doubly heavy mesons. Thereupon, perturbative corrections 
of order a s to the production amplitude were calculated in Ref. [l5|| increasing the cross 
section by a factor 1.8. On account of different values of relativistic corrections obtained in 
Refs.@,[il[ll[l3 and the importance of a relativistic consideration of the process e + +e 



J/ty + 1] C in solving the doubly heavy meson production problem, we have performed a new 
investigation of relativistic and bound state effects. As in our recent papers [3, 12 1, this 
investigation is based on the relativistic quark model which provides the solution in many 
tasks of heavy quark physics. In the above quoted papers Refs-QJ 1 171 1 we have demonstrated 
how the original amplitude, describing the physical process, must be transformed in order to 
preserve the relativistic plus bound state corrections connected with the one-particle wave 
functions and the wave function of a two-particle bound state. In the present paper we 
shall extend the method of Refs. 1(| 17 1 to the case of the production of a pair (V + V) of 
doubly heavy mesons containing quarks of different flavours b and c. In particular, we will 
consider the internal motion of heavy quarks in both produced pseudoscalar V and vector V 
mesons. The paper is organized as follows: In Sec. II we present the general formalism and 
the basic relations of our method which are required in order to formulate the relativistic 
amplitude for the production of doubly heavy mesons. In Sec. Ill we derive the analytical 
expressions for the corresponding cross sections and make numerical estimations exploiting 
the relativistic quark model. Conclusions and discussion of the results are given in Sec. IV. 
The construction of the quasipotential heavy quark distribution amplitude in the light-front 
variables is included in the Appendix. 



II. GENERAL FORMALISM 



The production of a pair of doubly heavy mesons in a color singlet model from e + e~ 
annihilation contains two stages after the transition of the virtual photon 7* into a quark- 
antiquark pair (Q1Q1). In the first stage examined on the basis of perturbative quantum 
chromodynamics (QCD), one of the heavy quarks (Qi or Qi) emits a gluon with sufficiently 
large energy ~ which then transforms to another heavy quark-antiquark pair (C^^)- 
These four quarks can combine with a definite probability into a pair of 5* - wave pseudoscalar 
(QiQ2)s=o an d vector (QiQ2)s=i doubly heavy mesons. The second nonperturbative stage 
of this process involves the formation of heavy quark bound states from heavy quarks. In 
the quasipotential approach to the relativistic quark model we can express the invariant 
transition amplitude for the described process as a simple convolution of a perturbative pro- 
duction amplitude of four heavy quarks T(pi,p 2 ; Qi, Q2), projected onto the positive energy 
states, and the quasipotential wave functions of a vector meson ^y{p, P) and a pseudoscalar 
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meson ^ v (q,Q) [18l ]: 

M(p-,p+,P,Q) = J j^^v(p,P) J ^^v(q,Q)-Mp 1 ,p 2 ;quq 2 )v(p + ) 1 a u(p.)^^ J 

(2) 

where p_, p + are four momenta of the electron and positron, r 2 = s = (p_ + p+) 2 , pi, P2 
are four-momenta of Q\ and Q2 quarks forming the vector doubly heavy meson V; q\, q 2 
are four momenta of Qi and Q 2 forming the pseudoscalar doubly heavy meson V. They are 
defined in terms of total momenta P{Q) and relative momenta p(q) as follows: 

E i2 Ml-m 2 21 +ml 2 
p 1 , 2 = r ]l>2 P±p, (p.P) = 0, ^2 = j^ = ^ -, (3) 

qi,2 = P x,2Q ± q, (q-Q)=0, p lj2 = ^- = ^ >-, (4) 

where My = m 1 + m 2 + Wy, Mp = m x + m 2 + IT-p are the masses of vector and pseudoscalar 
mesons consisting of heavy quarks. 

Different color-spin nonperturbative factors entering the amplitude A4(p-,p+, P, Q) con- 
trol the production of doubly heavy mesons. In this process the gluon virtuality is large 
k 2 ^> Aq CD , and the strong interaction constant is small a s C 1. There exist four produc- 
tion amplitudes in the leading order over a s as explained below in Fig.l. In a color singlet 
model the first amplitude P\^{p\,p 2 \ qi, q 2 ) takes the form: 

1x,pipi,P2;qi,q 2 ) = 16nas u 1 (p 1 )7 A , 7 ^ J 1 + m 2 1 ] . IpViiq^u^q^v^p^D^k). (5) 

3 (r — qi) z — mf + le 

The color factor 5uSkj(T a )ij(T a )ki/3 = 4/3 is already extracted in the amplitude (5). To 
calculate relativistic effects we have to keep all factors in the amplitude (2) with the relative 
motion momenta p and q. We have to take into account also the bound state corrections 
which are determined by the binding energies Wy and W-p. 




FIG. 1: The Feynman diagrams for the production of a pair of doubly heavy mesons (V + V) in 
e + e~ annihilation. The wave line corresponds to the photon and the dashed line corresponds to 
the gluon. Two other diagrams can be obtained by permutations. 

Clearly, in order to calculate the matrix element (2), one needs explicit expressions for 
bound state wave functions. We construct such wave functions by solving a relativistic 
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quasipotential equation to the desired accuracy in the center-of-mass (CM) frame. But in 
the matrix element (2) the final meson states have different total momenta P and Q. So, 
it is necessary to know how to transform the CM wave function to an arbitrary reference 
frame. The transformation law of the bound state wave function from the rest frame to the 
moving one with four-momentum P was derived in the Bethe-Salpeter approach in Ref. 
and in the quasipotential method in Ref. [2^. We use the last one and write the necessary 
transformation as follows: 

9ft (p) = D\' 2 > pa {RZ)D\ 12 ' W/3 «X(P), (6) 

*p (p) = % T {v)Dt 1A e \Rf P )Dt 1/2 > Tff « p ), 

where R w is the Wigner rotation, Lp is the Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix D 1/,2 (i?) is defined by 



'1 N 
,0 1, 



\D^(RW p ) = S- 1 (p h2 )S(P)S(p), (7) 
where the explicit form for the Lorentz transformation matrix of the four-spinor is 



S(P) = \ e(P)+m f 1 + / aP) ) ■ (8) 
V 2m V e(p) + mj y ' 

For further transformation of the amplitude (2) the following relations are applied: 

Stf(A)uj(p)= £ uZ(Ap)D%{R%), (9) 

(T=±l/2 

uJ(p)s£(A) = E Dt 1/2 (Rl P )Km- 

(T=±l/2 

Using also the transformation property of the Dirac bispinors to the rest frame 

Si(p)=Si(0) , === = , Pi = (ci,p), 
'2ei(p (ei(p)+mi) 



^(~P)= , = ^ M O), p' 2 = (e 2 ,-p), (10) 
/2e 2 (p)(e 2 (p) +m 2 ) 



we can introduce the projection operators ff p ' v onto the states (Q1Q2) in the meson with 
total spin and 1 as follows: 

ft* v = h(o)^(o)] 5=0il = 75(6*)^. (11) 

As a result the doubly heavy meson production amplitude from e + e~ annihilation takes the 
form: 

. , , ^ ^, 8n 2 aa s , . fl . , f dp ^o(p) 

Mi(p-,p+, P, Q) = -—^v{p + )^u{pJ) / \ 0Ky) x 

3s J ( 27r ) 3 v /2e 1 (p)(e 1 (p) + m 1 )2e 2 (p)(e 2 (p) + m 2 ) 

(12) 



5 



x 



dq 



^o(q) 



( 2?r ) 3 J2e 1 (q)(e 1 (q)+m 1 )2e 2 (q)(e 2 (q)+m 2 ) 



Sp{(p2 - m 2 )e*(l + vi)(pi + mi)x 



*7/ v ^ ^ 2 +m 2 \ ■ 7/?(gi - m i)7s(l + #2X92 + m 2 )7 i ,}-D /li ,(A;), 
(r — qi) 2 — mf + le > 

where the four- vectors e, p ij2 , q 1}2 are given by: 

e = L P (0,e)=(ev,e+^), 



(13) 



p li2 = S(Lp)ft.S-\L P ), S(L P )(1 ±-r°)S- 1 (L P ) = ±(t)i ± 1), t)i 



?i,2 = SiLcfifaS-^LQ), S(L Q )(1 ± 7°)^ 1 (^q) = ±(«2 ± 1), ^2 



P 

My' 

Mp' 



In order to make optimum use of the expression (12), we rearrange the bispinor contractions 
in the numerator of Eq.(12) and thus extracting in a more evident form the relative momenta 
p and q of heavy quarks: 



Mi(p-,p+,P,Q) 



x 
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dq 



%n 2 a.a. t 
~ 3s 

*?(q) 



v(p+)^u(p-) 



dp 



*£(p) 



( 27 r)3 y 



<siQ) (gl(p)+™-l) £2 (p) (£2(p)+m 2 ) 



(14) 



2mi 



(27T) 



^2 fa) (e2(g)+m 2 ) 
mi 2mi m 2 2m2 



f)l-l 



+ «1 
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X 
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+ Vl 



+ 
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2m\{e\{p) + mi) 2ui\ 



'V, 



2m 2 (e 2 (p) + m 2 ) 2m 2 
t — qi + mi 



I (l+i)i)x 



r — gi) 2 — mf + ie 



x 



? ! 2 
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5 



2mi(ei(g) + mi) 2mi 



75(1+^' 



«2 + 1 



l>2 



q^ 



q 



2m 2 (e 2 (q) + m 2 ) 2m 2 
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The production amplitude (14) contains relativistic corrections of two species. The first 
type of corrections are held in the quark interaction operator. They can be accounted for by 
means of the numerical solution of the Schrodinger equation with the relevant potential. The 
second type corrections are determined by several functions depending on the momenta of 
relative motion of quarks p and q including the gluon propagator D^ u (k), relativistic bispinor 
factors and the heavy quark propagator. At least, there exist bound state corrections in 
Eq.(14) which are related to Eqs.(3) and (4). 

Notice that there are several energy scales which can characterize the heavy quarkonium: 
the hard momentum scale ttiq (the mass of heavy quark), the soft momentum scale rriQVQ 
(vq is the heavy quark velocity in the bound state) and the ultrasoft momentum scale 
m Q v Q- We assume that the heavy quarkonium is a nonrelativistic system. This means that 
the following inequalities occur: uiq ^> rriQVQ ^> rriQVQ, ijiq ^> Aqcd, which we can exploit 
in the further study of the production amplitude (14). The expansion of basic factors over 
relative momenta p and q up to terms of the second order has the following form: 



1 

k 2 



1 



1 _ p 2 + g 2 2{p 2 (rp) + n 2 {rqj\ A[p 2 (r p) 2 



ri(rqf 



Zi 



Zi 



zl 



(15) 



6 



(r — q\) 2 — m\ 



1 

Z~2 



q 2 2(rq) 4(rg) 2 
Z 2 Z 2 



Z\ = i] 2 p 2 r 2 



Z 2 = r 2 -2 Pl (rQ)+piM£-mi, r = P + Q, 

where we have chosen the center-of-mass frame with the condition P + Q = 0; r = y/s is the 
total energy in the electron-positron annihilation. In the leading order, when one neglects 
relativistic and bound state corrections, one obtains Z\ = (1 — k) 2 s, Z2 = (1 — n)s, k = 
mi /Mo = mij (mi + m 2 ). So, the expansions (15) are well defined. Substituting (15) into 
Eq.(14) we have to combine the factors of the second degree over p and q. After averaging 
the obtained expression over the angle variables with the account of orthogonality conditions 
(3) and (4) 0, we can present the necessary correction in the integral form / dpp 2 ^ (p). 
Moreover, we performed the expansion of the amplitude (14) at p = q = over the binding 
energies W-p and Wy in the linear approximation in such a way that relativistic and bound 
state contributions can exist separately as follows: 



(T P ~# A 

xe ap \pv°vZe 



Mi(p-,p+,P,Q) 
M 



lTT 2 aa s 



v(p + )Yu(p-) 



(2tt) 3 
M 2 



(27T) 



3 <(q) x (16) 



sHl - k) 3 



2W V 
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M (1-k) v s 
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72M 2 k 2 (1 - k) 2 



-9 + 8k 4 + 



Mi 
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^ (1 _ 14 M + 10 M' 

s s 2 



72M n 2 K 2 (l - k) 2 



x 
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-9 + 16k 1 + —± - 16k 3 ^ 
Vs s 



' g _ 4Mo + UM[ 



2 M 



Mo 2M 2 



32k" — = -1 + — = - 



+ 



2Mf\ 
s 3 / 2 J 



Similar expressions occur for the other amplitudes contributing to the production process: 



x 



M 2 (p-,p+,P,Q) 
Mo L 



l7i 2 aa< 



s 3 (l - k) 3 
^2 



2W V M 2 . 
M 1 s ■ 



v( P+ )^u(p.) J t|^(p) / ^5*f (Q) >< M 
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M 2 
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Mqk(1 — k) s 
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dp 
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MnK 



■2 + 



2M 







72M 2 k 2 (1 - k) 1 



Mn 



Mn 



7—16—^—16- 



-160 : 



M 4 



Mn 



+4k 5 + 10-£ + 36 



s 

M 2 



2k + K I + 







Mn 



,3/2 



M 4 \ 

16^7 + 136— f + 



7 

e+ + e~ -> J/* + ?? c 




7 8 9 10 11 

e+ + e" B+ + B*" 




e+ + e _ — > T + rj h 






008 









006 









004 


/ / 
/ / 
/ 







002 


/ 

/ 

/ 

/ 









/ 








20 21 


22 23 24 25 26 



s 



FIG. 2: The cross section in fb of e + e~ annihilation into a pair of doubly heavy mesons with 
opposite charge parity as a function of the center-of-mass energy yfs (solid line) . The dashed line 
shows the nonrelativistic result without bound state and relativistic corrections. 
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It is important to emphasize that we carried out formal expansions of the integrands in 
the amplitudes Aii (16)-(19) over quantities p 2 /m 2 2 , q 2 / m i2 assuming that the relative 
momenta of heavy quarks in the bound state are small as compared with their masses. This 
means that the formally divergent integrals / dpp 2 ^(p) and / rfqq 2, I , (q) have to be defined 
by a suitable regularization procedure. Their numerical values will be directly determined 
by the properties of the wave functions of the heavy quark bound states. 
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III. CROSS SECTION OF PSEUDOSCALAR AND VECTOR DOUBLY HEAVY 

MESONS PRODUCTION 



Using the relativistic amplitudes Aii obtained above, we can calculate the production 
cross section for doubly heavy mesons in e + e~ annihilation as a function of the center-of- 
mass energy y/s. For this purpose, the sum of amplitudes Aii (16)-(19) can be represented 
as follows: 

16 2 M y/4M V M V 

— 7T n 

27 



V + V) = Z V :Z ^(0)^(0) x 



M(e + e- 
xu(p+)/u(p_)e ffpA/ 9<u£e* A 



S 6 K 5 (1 — k) E 

K 3 Qia a2 7i + (1 - nfQ 2 a sl T 2 



(20) 



where ot s i = a s (4m 2 ), a s2 = a s {Am'^), <5i,2 are the electric charges of heavy quarks, ^'^(O) 
are the wave functions for the relative motion of heavy quarks in the vector and pseudoscalar 
meson at the origin in its rest frame, 
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Ma ' M 

where (p 2 ), (q 2 ) are quantities determining the numerical values of relativistic effects con- 
nected with the internal motion of the heavy quarks in vector and pseudoscalar doubly 
heavy mesons. Note that they are not equal to the matrix elements / rfp^Q (p)p 2 \E f g (p) and 
/ dq 1 4 r o'(q)q 2 ^o'(q) and are discussed at the end of Sec. III. 

Performing standard algebraic calculations with the squared modulus of the amplitude 
\Ai\ 2 by the use of the system Form j^jj, we obtain the following differential cross section: 



do 



An 3 a 2 M 2 \^ 



W(o)p 



dcosfl 729M v M p s 8 k 10 (1 - k) 10 

f [s 2 - {M v + M v ) 2 } [s 2 - (M v - M v f 



K 3 Qia s2 T\ + (1 - K) 3 Q 2 a sl T 2 

3/2 



(23) 



1 + cos 2 fl , 



where 9 is the angle between the momenta of the incoming lepton and the outgoing doubly 
heavy meson. The total cross section for the exclusive production of pseudoscalar and vector 
doubly heavy mesons in e + e~ annihilation is then given by the following expression: 



a is 



327r 3 a 2 M 2 |tf£(0)| 2 |tf£(0)| i 
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(24) 



1 - 
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v, 



1 - 



(My - M 
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21 1 3/2 



If we neglect relativistic and bound state effects in Eq.(24), that is (p 2 ) = (q 2 ) =0, Mp = 
My = M , then we obtain unr(s) coinciding with the analytical result of Ref.|22j. To 
estimate how the relativistic and bound state corrections can change the nonrelativistic 
cross section <Jnr(s), we have used definite numerical values of a number of parameters 
entering in Eqs. (21)-(24). They can be determined on the basis of the relativistic quark 
model 2^, 24], 25| as we shall discuss in the following. 
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FIG. 3: The cross section in fb of e + e~ annihilation into a pair of S'-wave doubly charm heavy 
mesons with opposite charge parity as a function of the center-of-mass energy s (solid line). The 
dashed line shows the nonrelativistic result without bound state and relativistic corrections. 
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In the quasipotential approach the bound states of heavy quarks are described by the 
Schrodinger type equation 23] 

(b\M) p 2 \ T< , . f dq 



\ 2/2 R 2jj, R J J (2tt 

where the relativistic reduced mass is 



*o(p) = / ?A^(P> q> M)* (q), (25) 



ExE 2 M 4 -(m 2 -m 2 ) 2 
m ~ E^TW 2 ~ ' (26) 

and the particle energies E\, E 2 are given by 

= M 2 - m 2 2 + m\ = M 2 -m\ + m\ 

1 2M ' 2 2M ' [ ' 

Here M = Ei + E 2 is the bound state mass, are the masses of heavy quarks (Qi and 
Q2) which form the meson, and p is their relative momentum. In the center of mass system 
the relative momentum squared on the energy surface M = E\ + E 2 reads 

2 [M 2 -{m x + m 2 f\[M 2 -(m x -m 2 f\ 

b (M) " UP • (28) 

For small binding energies it is approximately equal to 2^rW. Notice that the kernel 
V(p, q, M) in Eq. (j25|) is the quasipotential operator of the quark- ant iquark interaction. 
The construction of the quark interaction operator is discussed permanently during the last 
decades jll 0, 0, . A new stage in solving this problem came when nonrelativistic field 
theories were introduced for the study of nonrelativistic bound states 27, 2^, 2^. Within an 



effective field theory (NRQCD) the quark-antiquark potential was constructed in Refs.|2i 
I30I l3~D ] by the perturbation theory improved renormalization group resummation of large 
logarithms. On the other hand, in the quasipotential quark model the kernel V(p, q, M) 
is constructed phenomenologically with the help of the off-mass-shell scattering amplitude, 
projected onto the positive energy states. The heavy quark-antiquark potential with the 
account of retardation effects and the one-loop radiative corrections can then be presented in 
the form of a sum of spin-independent and spin-dependent parts. The explicit expression for 
it is given in Refs.jiil. l25l|. Taking into account the accuracy of the calculation of relativistic 
corrections to the cross section (24), we can use for the description of the bound system 
(Q1Q2) the following simplified interaction operator in the coordinate representation: 

V(r) = ~^^- + Ar + B, (29) 
where the parameters of the linear potential A = 0.18 GeV 2 , B = —0.3 GeV, 

47T 2 

asW = Wp)' A, = 11 -3"/- (30) 
Here nt = 3 is the number of flavours and a = 2mir " 2 is a renormalization scale and 

J r" (7711+7712) 

A = 0.168 GeV. All the parameters of the model like quark masses, parameters of the linear 
confining potential A and B, mixing coefficient e and anomalous chromomagnetic quark 
moment k entering in the quasipotential V(p, q, M) were fixed from the analysis of heavy 
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quarkonium masses and radiative decays Solving the Schrodinger- 

like quasipotential equation with the operator (29) we obtain an initial approximation for 
the bound state wave functions in the case of (cc), (be) and (bb) systems. For numerical 
estimations of relativistic and bound state effects in the production of heavy mesons in the 
e + e~ annihilation we need the values of the wave functions at the origin, the bound state 
energy and the parameter of relativistic effects: / p 2 \E r ^' v (p)dp/ (27r) 3 . This integral diverges 
at high momentum. Using the Schrodinger-like equation (25), we can express it through the 
value of the potential V(r) at the origin. Different regularizations for it are discussed and 
already applied to a number of tasks concerning the production and decay of bound states 
^"2, 33, 3J, 35[. A recent discussion of the method for the calculation of an order-f 2 matrix 



element is presented in Ref. |2a]. To estimate the numerical value of relativistic corrections 
in the production processes, we follow the prescription of dimensional regularization, where 
the scaleless momentum integral JV(p — q)\E r (q) -Jjp related to the problem vanishes 
Then we can express the necessary quantity in the form: 

P 2 ) D « s W) I (SpP 2 *°<P> = + 2,n\B\. (31) 

Let us note, that in the case of the Coulomb interaction in QED the dependence of the 
relativistic parameter (31) on the principal quantum number n is known in the analytical 
form F° r the quark bound states the dependence on the quantum number n is deter- 
mined numerically. Solving the quasipotential equation (25) with the potential (29) [38j|, we 
obtain the energy spectrum W of the heavy quark system and the numerical values of the 
parameter (31) for the bound states (cc), (bb) and (6c) which are presented in Table I. Heavy 
quark symmetry predicts that the wave functions of the vector and pseudoscalar states are 
different due to corrections of order v%. The analogous statement is valid for the parameter 
(p 2 ). Neglecting in this study the corrections of order 0(vq) in the production rates we 
write in Table I equal values of (p 2 ) for V- and P-mesons. The theoretical uncertainty 
of the obtained value of the parameter (p 2 ) in Table I is determined by perturbative and 
nonperturbative corrections to the quasipotential 0, 24]| and is not exceeding 30%. We 



presented in Table I the estimation of the pseudoscalar and vector wave functions at the 
origin including effects of order 0(vq). For this aim first of all, we obtained the central equal 
values of \I/y(0) and \l/-p(0) solving the Schrodinger-like equation (25). Then we took into 
account the spin- dependent corrections of order 0(vq) to the wave functions at the origin 
which were calculated in the framework of nonrelativistic QCD (see 



IV. SUMMARY AND DISCUSSION 

We have investigated the relativistic and bound state effects in the production of doubly 
heavy mesons on the basis of the perturbative QCD and the relativistic quasipotential quark 
model. Using the factorization hypothesis we keep systematically all relativistic corrections 
of the second order in the relative velocity of heavy quarks and bound state effects in the 
production of the S'-wave pseudoscalar and vector mesons from e + e~ annihilation. 

Let us summarize the basic peculiarities of the calculation performed above. 

1. We obtain the cross sections for the production of a pair of S'-wave doubly heavy 
mesons with opposite charge parity containing b and c quarks from e + e~ annihilation. 
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TABLE I: Basic parameters of the relativistic quark model 



State 


Particle 


Mass, GeV 


a s 


Bound state 


tf(0), GeV 3 / 2 




n 2S+1 Lj 




[24, 42] 




energy, GeV 




GeV 2 




Vc 


2.980 


0.314 


-0.120 


0.28 


0.6 


l 3 Si 




3.097 


0.314 


-0.003 


0.26 


0.6 


2 x 5o 


rfc 


3.594 


0.314 


0.494 


0.25 


1.4 






3.686 


0.314 


0.586 


0.23 


1.4 




Bt 


6.270 


0.265 


-0.160 


0.36 


0.7 


13 C 


D c 


6 332 


265 


-0 0Q8 

U.V7 JU 


34 


7 


2 1 S 




6.835 


0.265 


0.405 


0.34 


1.8 


2 3 5i 


b' c *- 


6.881 


0.265 


0.451 


0.32 


1.8 


l x 5 


Vb 


9.400 


0.207 


-0.360 


0.55 


0.9 


l 3 5i 


T 


9.460 


0.207 


-0.300 


0.53 


0.9 


2^0 


v' b 


9.993 


0.207 


0.233 


0.45 


2.9 


2 3 5i 




10.023 


0.207 


0.263 


0.43 


2.9 



TABLE II: Comparison of theoretical predictions with experimental data. 



State 
H1H2 


CBaBar x 
BrH 2 ^charged>2 

m m 


CBeiieX 
BrH 2 -*charged>2 

(fb) [6] 


[14] 


&NRQCD 

(fb) [8J 


[10] 


_ 


Our result 
(fb) 


tt(lS)7fc(lS) 


17.6 ± 2.8±|:f 


25.6 ±2.8 ±3.4 


26.7 


3.78 


5.5 


7.4 


7.8 


*(2S) Vc (lS) 




16.3 ±4.6 ±3.9 


16.3 


1.57 


3.7 


6.1 


6.7 


*(15)t/ c (25) 


16.4 ± 3.7±|;^ 


16.5 ±3. ± 2.4 


26.6 


1.57 


3.7 


7.6 


7.0 


^>(2S) Vc (2S) 




16.0 ±5.1 ±3.8 


14.5 


0.65 


2.5 


5.3 


5.4 



2. All possible sources of relativistic corrections including the transformation factors for 
the two quark bound state wave function have been taken into account. 

3. We have investigated the role of relativistic and bound state effects in the total 
production cross sections using predictions of the relativistic quark model for a number of 
parameters entering in the obtained analytical expressions. 

The results of our calculation of the cross section (24) presented in Figs. 2-3 evidently show 
that only the relativistic analysis of the production processes can give reliable theoretical 
predictions for the comparison with the experimental data. It follows from Fig. 3 that with 
the growth of the quantum number n the nonrelativistic approximation doesn't work near 
the production threshold because the omitted terms in this case have the same order of the 
magnitude as the basic terms. 

As we have already mentioned, the experimental results for the production of J/\I/ + r\ c 
mesons measured at the Belle and BaBar experiments differ from theoretical calculations 
in the framework of NRQCD. The experimental data on the production cross sections of 
a pair of S-wave charm mesons are presented in Table II. The numerical value for the 
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cross section of J/*f> + r\ c production at sfs = 10.6 GeV, obtained on the basis of Eq.(24) 
amounts to the value 7.8 fb without the inclusion of QED effects. In this case relativistic 
and bound state corrections increase our nonrelativistic result by a factor 2.2 (cf. dashed 
lines in Figs.2,3). Accounting slightly different values of several parameters used in our 
model in the comparison with Ref.jSj (the mass of c quark, the binding energies W-py), 
we find agreement of our calculations with the results of Ref.jsf for the production of the 
charmonium states, if relativistic corrections are taken into account (see the seventh column 
of Table II) . Taking in mind also the calculation GH , which includes additional perturbative 
corrections of order a s , we observe the convergence between the experimental data and 
theoretical results obtained on the basis of approaches combining nonrelativistic QCD and 
the relativistic quark model. Our results should be useful for the future comparison with the 
(bb) and (cb) meson production measurements. Similarily, we can examine the production 
rates for doubly heavy mesons including the P- and .D-wave states. The experimental data 
for the production of the mesons Xco + J/^ are obtained already in Refs.0, 0|. The work in 
this direction is in progress. 

There is an essential difference in the numerical results obtained in this work and in 
Refs. HI QUI on the basis of the light-cone formalism. Our results lead to the increase 



of the cross section for the production J/^f + i] c only by a factor 2 ~ 2.5 in the range of 
center-of-mass energies yfs — 6 4- 12 GeV but not to an order of magnitude at y/s = 10.6 
Gev as in 1 141] . While the cross section (24) depends on the choice of the c-quark 

mass, strong coupling constant a s , the meson wave functions at the origin, binding energies 
W-py and the relativistic parameter (31), a possible growth of the theoretical value (24) 
doesn't solve the problem. So, there exist at least two questions which could be discussed 
regarding the performed calculation. 

The first question refers to the comparison of our calculation with the light-cone approach 
to the same problem in Refs. [l2l Il3l Since the main part of the investigation was 

connected with the relative motion of heavy quarks in the meson, it is necessary to find 
the wave function ip(p±,x) which describes bound states in such an approach and satisfies 
the quasipotential wave equation in the light-front formalism ji^, 44[ ^ (see Appendix A). 



This equation gives the wave function of a bound state in an arbitrary Lorentz reference 
frame. The transformation property for the wave function from the arbitrary frame to the 
frame, in which the total transverse momentum of the two-particle bound state is equal to 
zero is determined by Eq.(A3). But contrary to our approach to the production of doubly 
heavy mesons, accounting the relative motion of the heavy quarks, the quark transverse 
momentum inside the heavy quarkonium was neglected in Refs. [12I H3I 14l . So, the basic 
difference between our method and the light-cone formalism consists in the fact that in 
Ref.0 E! the effects of the wave function ip(x,p±) transformation from the frame 



P_i_ = to the moving one with the transverse momentum Pj_ 7^ were not taken into 
account. They lead to terms proportional to the meson transverse momentum squared Pj_ 
in the production amplitude and might correct the total result. 

Finally, another important question regards the accuracy of the calculation performed in 
this paper and the numerical estimation of the next term in the expansion of the produc- 
tion amplitude over the relative velocity of the fourth order. Recently, corrections of order 
0(f 4 ) were studied in the S— and P-wave quarkonium decays in Refs. 4(| 47]. There it was 
demonstrated that the v expansion converges well for the decays of heavy quarkonium. Note 
that a numerical estimation of the parameter (p 4 ) can be obtained as in Eq.(31) by using 

4 / m 4 \ _ {W+\B\f 



the dimensional regularization technique |35| : \P / m Q/ DR = — m L ■ When we consider 
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the production of IS*- wave charmonium states with the binding energy W -C rriQ, then the 
numerical values of the parameters (p 2 ) DR , (p 4 )dr are sman - I n the production of excited 
S-wave states the binding energy increases what leads to the growth of (p 2 ) DR , (p 4 )dr- 
it follows from the results of Table I that (v 2 ) Dij , m 0.6 for the 2S'-charmonium states and 
(v 2 ) M fa 0.7 for the 2S (be) meson. Obviously, the convergence of the expansion becomes 
worse in this case. The calculation of the coefficient in the term (v 4 ) would deserve addi- 
tional investigation. On the whole, we estimate the theoretical uncertainty of the obtained 
production rates due to inaccuracies in the determination of the meson wave functions at 
the origin and the omitted terms in the expansion of the relativistic amplitudes as the 40% 
for the 15-states and the 70% for the 25-states. 
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APPENDIX A: QUASIPOTENTIAL LIGHT-CONE DISTRIBUTION 

AMPLITUDE 

In this Appendix we consider the construction of the heavy quark light-cone distribution 
amplitude on the basis of the quasipotential equation derived in light-cone variables. 

The quasipotential wave equation for the two-body bound state in the light-front formal- 
ism takes the form ji^, 44, 45[: 



P 2_ 



(p±+ (|-x)P ± ) 2 + m2 (p ± + (|-x)P ± ) 2 + 



rrin 



X 



X 



1pp(P.L,x) 



(Al) 



dy 



V(P, x, p ± ; y, q±)ip P {q±, y), 



lo 2/(1 - y) J 2(2tt)3 

where p (p±) and P (Pj_) are total (transverse) momenta of relative motion of quarks and 
the meson. The variable x is introduced in the following way: 



x = - + p± = Po ± pa, P± = Po± P3- 



(A2) 



This equation allows to find the two-particle wave function in an arbitrary Lorentz reference 
frame. The transformation of the wave function from the arbitrary frame to the frame, in 
which the total transverse momentum P^ = is the following: 



ip P (x, p ± ) = ij F± =o(x, p ± + (1/2 - x)P_l) 



(A3) 
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The case of spin particles is discussed in Refs.j45, 48[. In the equal mass case (mi = m 2 = m) 
in the frame where Pi = we get: 



2 p ± z + nr 



M - 



x(l — x) 



ip{p±,x) 



dy 



o y(l-y)J 2(2vr)3 



:V(x,p ± ;y,q±)ip(q±,y)- (A4) 



2.5 



1.5 



0.5 




FIG. 4: The light-front distribution amplitude over the momentum fraction x (solid line). The 
dashed line shows the phenomenological function from Ref. }13j]. The thin solid line shows the 
result obtained in our model by means of the relations (44) , (49) of Ref . j4jj with the c-quark mass 
m c = 1.55 Gev. 

In order to discuss the realistic wave function ip(p±,x) there is the need to transform 
Eq.(A4) to the ordinary form (25), because we know the potential (29) for it which describes 
the energy spectrum with sufficiently high accuracy. Let us introduce the relative momentum 
p = (p_i_,£>3) = (psin^cos^psin^cos^pcos^): 



x — 



P3 



1 



- y 'p± 2 + m 2 , x = -[l 

Jx(l — X) z \ 



P3 



so that 



dxd 2 p± 



e(p) 
2x(l — x 



1 / pcos6 \ 

2 l 1+ V^ + mV ' 



dp. 



(A5) 



(A6) 



^Jp 2 + m 2 

As a result, we obtain from Eq.(A4) the ordinary quasipotential equation (25), which con 
tains the wave function ^o(p) connected with the initial one by the following relation: 



*o(p) 



4x(l — x) 



nl/4 



p± 2 + m 2 



ip(px,x) 



(A7) 



The ordinary normalization condition for the wave function ^o(p) can De rewritten for the 
function ip(p±,x) as follows: 



dx 

o x(l — x) 



dp± 



[P±,x) 



(A8) 
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Then we determine the light-cone distribution amplitude by means of the quasipotential 
light-front function ip(p^,x): 

= M J ^( P± , x)d 2 P± = M • 8tt • x(l - x) j™ r^j^ jp 2 + m 2 ) 1 ^ (p)dp. (A9) 

The integral function in Eq.(A9) differs by the power of the relativistic energy e(p) (e 1//4 (p) — > 
e 1 ^ 2 (p)) and the function x(l — x) from the expression in Ref.[49j. The normalization factor 
M is determined by the condition: Jq 1 4>(x)dx = 1. The plot of the function 4>{x) and the 
comparison with the corresponding amplitudes in Refs. (la. |49| are presented in Fig.4. All 
these functions, obtained in the frame where the meson transverse momentum = 0, have 
the similar shape. But our function (solid line in Fig.4) has a larger maximum and falls 
quicker to the end points x — 0, 1. To derive the light-cone distribution in the arbitrary 
reference frame we can use Eqs.(A3), (A9). The shift of the transverse momentum p± — > 
p_i_ — Q — x^j Pj_ shows that the quark light-cone distribution doesn't change the shape in 
the arbitrary frame with 7^ 0. So, the light-cone amplitude transformation from dashed 
line of Ref.hJ to our solid line in Fig.4 decreases the production rates of the charmonium 
states in the light-front approach. 
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